Conductance fluctuations as a tool for investigating the quantum modes 
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Recently it has been observed that the conductance fluctuations of atomic size gold contacts are 
suppressed when the conductance is equal to an integer multiple of the conductance quantum. The 
fact that these contacts tend to consist exclusively of fully open or closed modes has been argued 
to be the origin for this suppression. Here, the experiments have been extended to a wide range 
of metallic elements with different chemical valence and they provide new information about the 
relation between the mode composition and statistically preferred conductance values observed in 
conductance histograms. 
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INTRODUCTION 



Manipulation and characterization of atoms and 
atomic size metallic constrictions has recently, become 
available through the development of the STMId. An al- 
ternative tool, for creating stable and clean atomic size 
metallic contacts, js the mechanically controllable break 
junction (MCB)aa. For the characterization of these 



systems, measurements of the electrical conductance are 
widely employed. This, as a result of the ease with which 
they usually can be obtained. The framework within 
which one should describe the conductance of such small 
contacts, which in the case of metals have dimensions of 
order of tbf; Fermi wavelength, is the Landauer-Biittiker 
formalismQ. In this formalism the conductance in the 
contact is described by N channels, determined by the 
narrowest cross-section of the constriction and the Fermi 
wavelength. Each channel has a transmission probabil- 
ity T„ with a value between and 1. The total conduc- 
tance is given by G = {2e^/h) '^n- For an adiabatic 
constriction in a free electron gas, the conductance in- 
creases stepwise with quantum units of the conductance 
(Go = 2e^//i) as the channels operL.one by one while 
increasing the constriction diameteitl. However, when 
one pulls apart a metallic atomic size contact, neither 
the diameter nor the conductance of the constriction de- 
creases smoothly. Instead, a series of steps (of order Gq) 
and plateaus are observed in the conductance on elonga- 
tion of the contact. The sequence of steps and plateaus 
is different each time the contact is pulled apart. The 
steps correspond with atomic reconfigurations and the 
plateaus with elastic deformation of the contactfl. It is 
tempting, but in principle incorrect, to assume offhand 
that these conductance measurements on atomic necks 
simply probe a series of discrete diameters of a free elec- 
tron gas. It may work for some metals (we will show that 
for sodium this is nearly the case), but a correct general 
description of the conductance of metallic point contacts 
consisting of even in the simplest case, a singlepa|tpm, has 
to consider the chemical valence of this atomQ^El. Com- 
pare for instance a single atom contact of the monovalent 



s-metal gold with the trivalent sp- metal aluminum. In 
the former case the conductance is carried by a single 
channel with conductance close to Gq. A single atom 
of aluminum on the other hand, with three conduction 
electrons in the 3s and 3p shells, also has a total conduc- 
tance close to Go but allows three partially transmitting 
channels. 

Recentl};^ a new technique, making use of conduc- 
tance fluctuations, has bfi©¥ presented which does not re- 
quire superconductivitylsllj to obtain information about 
the conductance modes contributing to the conductance. 
First results on gold contacts with conductance up to 
5 Go have shown that for this s-metal, once a channel 
(with number n) is partially open it tends to fully open 
before a next (n-l-1)*'' channel starts to contribute signifi- 
cantly. This interpretation was confirmed by an indepen- 
dent technique which consisijs^af measuring the shot noise 
in the point contact currcntO. In this paper we present 
a more complete argumentation of the theory together 
with new measurements of the conductance fluctuations 
on copper, silver, sodium, aluminum, niobium and iron 
and discuss what information on the channel transmis- 
sions can be extracted from our results. We will also 
compare our conductance fluctuation results to recently 
published measurements on the thermopower of atomic 
size contacts and show that both measurements can be 
related without any free parameters. 

Fluctuations in the conductance with bias voltage haiie 
previously been observed in larger ballistic contactst3 
and have an origin analogous to the universal conduc- 
tance fluctuations (UCF) measured in diffusive wiresli^. 
The interesting new aspect of such fluctuations in quan- 
tum point contacts is that their r.m.s. amplitude depends 
on the transmission probability of the channels contribut- 
ing to the conductance. The underlying principle of this 
effect can be understood by considering a contact with a 
single conducting mode having a finite transmission prob- 
ability T, described by transmission and refiection coef- 
ficients t, t' , r and r' (coming in from left and right, re- 
spectively), with = |i|2 = T and jr'p = |rp = 1 - T. 
As illustrated in Fig. , electron waves transmitted by the 
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contact with amplitude t, and scattered back towards the 
contact through diffusive paths in the bank with proba- 
bihty amphtude a, have a probability amplitude r to be 
reflected at the contact. This wave interferes with the 
directly transmitted partial wave and modifies the total 
conductance, depending on whether the resulting inter- 
ference is constructive or destructive. A similar contri- 
bution comes from the trajectories on the other side of 
the contact. The interference terms will be sensitive to 
changes in the phase accumulated along the trajectories, 
which is determined by the electron energy and the path 
length. The fluctuations in the conductance as a function 
of bias voltage are thus the result of the change of these 
phase factors by the increase in the kinetic energy of the 
electrons by an amount eV. What is immediately appar- 
ent from the principle illustrated in Fig. is that when the 
coefficient T is either 0, or 1, the interference and thus 
the amplitude of the fluctuations vanishes. This suppres- 
sion of conductance fluctuations at quantized values has 
been noted in numerical simulations of quaptum point 
contacts containing disorder by Maslov et al!t3 
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FIG. 1. Schematic diagram of the configuration used in 
the theoretical analysis. The dark lines with arrows show 
the paths, which interfere and contribute to the conductance 
fluctuations in lowest order. 



Each time when the contact is opened, and closed again 
to sufficiently large conductance values, random atomic 
reconfigurations take place, leading to a completely new 
set of scattering centers. The statistical results of many 
different contacts can hence be interpreted as the ensem- 
ble average over defect configurations. With this tech- 
nique we have studied the average properties of the con- 
ductance modes for different materials and their relation 
with the statistically preferred condiietaace values ob- 
served previously by various authorstaiHj through the 
measurement of so called conductance histograms. 

Typically when studying conductance fluctuations one 
measures the differential conductance over a wide range 
of bias voltage or magnetic field. Here, on the other 
hand, we measure the first and second derivative of the 
current with respect to voltage of atomic size contacts. 
The first derivative gives us the conductance, and the 
second derivative is a parametric derivative of the con- 



ductance. The latter can roughly be seen as a measure 
for the amplitude of the fluctuations with voltage. We 
measure these quantities for a large number of different 
contacts, thus determining effectively an average over an 
ensemble of scattering configurations. This measurement 
method is preferable as it is much faster than measuring 
the conductance as a function of bias voltage directly, 
and hence allows the experimental determination of the 
average properties of many contacts within a reasonable 
time scale. 



THEORY 

In this section we will give a more detailed description 
of the phcnomcnological theory presented in Rcf. |9|. In 
our model for a metallic constriction we divide the con- 
ductor into three separate regions: A region, small on the 
scale of the scattering lengths involved, centered around 
the narrowest part of the conductor, which we describe as 
ballistic. On either side of this ballistic region we consider 
a diffusive region, characterized by a mean free path, Ig. 
In order to make the geometry of the contact more realis- 
tic, we assume a conical shape for these diffusive regions, 
with opening angle 7 (see Fig. ). The probability ampli- 
tudes for scattering from any incoming mode on the left 
to any outgoing mode on the right side of this ballistic 
section (henceforth referred to as the 'bare contact') is 
described in terms of the matrices of transmission t, t' 
and reflection r, r' when coming from the left and right 
respectively. 

where M is the transfer matrix and v, Or, h, oi are the 
vectors of the incoming and outgoing waves for the right 
and left hand side respectively. The matrix of trans- 
mission probahilities is given by T = tt^, which can be 
diagonalizedEifEj. For a narrow constriction, most of the 
diagonal elements will be zero. The number of conduct- 
ing modes N and their transmission probability is given, 
by the non-zero diagonal elements, T„ = = 
n = l,2,...,iV. The reflection probability of mode n is 
given by, i?„ = |r^| = |r^f | = 1 — r„. In the simple free 
electron gas model, the number N is determined by the 
width of the narrowest part of the contactQ and by the 
Fermi wavelength. Note, however, that in principle it is 
not restricted to any particular model. The values of the 
r„'s are somewhat influenced by our arbitrary choice of 
the boundaries between the ballistic and diffusive regions. 
For this influence to be small, the distance L between the 
center of the contact and the boundaries should be large 
on the scale of the contact diameter. On the other hand, 
in order to be able to neglect fluctuations of r„ on the 
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scale of the applied voltage, V, we require L <C hvp/eV. 
For metallic contacts we can typically take L ~ 1 nm. 

The left and right banks are also described in terms of 
transfer matrices M; ,,7 similar in form to the one used 
for the bare contact in Eq.|l|. We define in this case 
and Sir as the transmission matrix and return amplitude 
matrix for the right bank and t; and a; the correspond- 
ing ones for the left bank. The elements of the return 
amplitude matrices a;^^,^ and ar„„ , which scatter a wave 
from mode m to mode n in the left or right hand side dif- 
fusive section, are expected to be small compared to 1. 
It is hence, a reasonable approximation to calculate the 
total transmission probability to first order in the return 
amplitude matrix elements only. The return amplitudes 
are energy dependent but this will not be made explicit 
until this dependence becomes relevant in Eq. ^ The to- 
tal transmission matrix for the two banks and ballistic 
constriction combined can be written as. 



tt = ((M,MM022) 



(2) 



t' 



Since the return amplitudes will usually be small we can 
set ti = tr — I, the identity matrix. With this assump- 
tion and the fact that we are only calculating to lowest 
order in ai^rmn neglecting corrections to the to- 

tal conductance and higher order contributions to the 
conductance fluctuations. These corrections will be dis- 
cussed at the end of this section. 

In order to calculate the conductance fluctuations, 
we will make an expansion of the total transmission 
probability to lowest order in ai^r„„- Using tit| = 
{{t^^y {t^^))~^ and assuming the matrices t, t', r and r' 
are already in diagonal form, the trace of ttt| to lowest 
order in a;^r„„ can be approximated by the sum of the 
inverse of the diagonal components of (t^^)'l'(t^"'^), 
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« E^"(l + 2Re(a,„„r; + r„o,„J). (3) 

Tl=l 

We assume that the Boltzmann constant times the tem- 
perature, kBd, is much smaller than the energy scale of 
the applied voltage, eV, (in accordance with the situation 
in our experiment) so that we can take the zero temper- 
ature approxiHiation. The current is then determined by 
Eq.| through^. 



2e2 
h 



/ = — / Tv[tttl]dE. 



(4) 



The fluctuations in the conductance are described by 
SG = G - (G), with G = dl/dV. Combining this with 
Eq.| and Eq.g and using (G) = i2e^ /h)J2^^iTn to 
first order in a/ , leads to, 



(5) 



c2Re[r„a;„„(eV" - E) + ar,„^{-E)r'„]dE 



= - E — ^"2Re[r>,„„(el/) - r„ai_(-eV)]. 

n— 1 

The return amplitudes contain random phase factors of 
the form exp[— ^(i?;^ ± eV)T/h], where r is the traversal 
time for a particular trajectory. Averaging Eq. || over the 
ensemble of defect configurations will give a zero result, 
((5G) = 0. For the correlation function of the conduc- 
tance as a function of voltage, however, we obtain a fi- 
nite contribution. In the product, only terms of the form 
a;^r„„(-E'i)a^^ {E2) have a chance to survive the averag- 
ing. In addition, diffusion in the left and right banks is 
uncorrelated, so that products of ai^^{Ei) and ar^^{E2) 
average to zero, 

N , -2. \ 2 

(5G(el^i)5G(ey2)) =E (x^'V ^" 

n— 1 ^ ^ 

X (2ReK„„(eyi)<_(ey2) + a/„„(-eyi)a;*_(-ey2)]) • 

At this point we assume that the average properties 
of the scattering on both sides of the contact, for all 
the mode indexes, are the same. Further, we pro- 
pose that (a/,r„„(£'i)a*^^^(£'2)) can be expressed as 
/o°° -fc!(''')e~*^^^~^^-''^^''dT, where r is the time required 
for the completion of a classical diffusive trajectory, and 
Pel if) is the classical probability distribution to return to 
the contact at this time. The classical return probability, 
Pc/(t), can be obtained by considering an electron being 
injected from the ballistic central section of the contact 
into the diffusive region at the left or right. When we 
take the interface between the ballistic and diffusive re- 
gions to be at a small distance L from the contact center, 
then after a given time r, with Dr 3> i^, the probability 
distribution to find the electron at a distance r > L is 
given by the classical result, 



p(t, r)c?T 



(1 -cos7)(47ri:)T)3/2 



where D = vple/'i is the diffusion constant with l,, the 
mean free path for elastic scattering. Here we have as- 
sumed that the diffusive region has the shape of a cone 
with opening angle 7 (Fig. ) and that only a small num- 
ber of channels are transmitting, so that most electrons 
entering the ballistic region are reflected. The probability 
per unit time to find the particle back in a disk of radius 
a and thickness dx at the entrance of the ballistic region 
is p(r, L)'Ka^dx. The average time the particle spends in 
this volume is dx/{vx) = V^dx/vp. The probability that 
the particle moves towards the contact instead of away 
from it is 1/2, and we assume that it has equal proba- 
bility to enter into any of the Nt, modes available at the 
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entrance of the ballistic section, where Nj, = {kpa/2)'^. 
Thus, the probability per unit time to return to the con- 
tact after a time r, into a given mode, is 



Pciir) 



2V3^fc2 (i:>T)3/2(l - COS7) ' 



(7) 



This distribution should further be multiplied with a fac- 
tor e^^'^'t' in order to take into account the probability that 
within a typical time the electron undergoes scattering 
which destroys phase memory. Combining these expres- 
sions we obtain, 



{SG{Vi)SG{V2)) 



n=0 ^ 



Tn) 



X / Pc/(t) cos 
/o 



(8) 



In our measurements we really measure {{dG/dV)'^) 
with a fixed modulation voltage rather than 
{SG{eVi)dG(eV2))- This can easily be corrected by dif- 
ferentiating Eq. with respect to Vi and V2 , and then 
setting V = Vi ^ ¥2- In the limit of V — > and using the 
above approximation for the average return probability 
we obtain. 
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/Z^f^kl 1 - C0S7 D3/2 



Evaluating the integral and taking the square root, re- 
sults in, 
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hkpVF^/^ — COS7 
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The typical time scale on which a first collision takes 
place is the elastic scattering time = le/vp- 

We have conducted our experiment by measuring the 
first and second derivative of the current with respect to 
voltage. The amplitude of the applied modulation volt- 
age was responsible for the energy cut-off rather than 
the dephasing time as we have assumed above {eVmod > 
h/r^). Using the derivation in the Appendix, which in- 
corporates this finite mpdulation voltage into the theory 
gives us the final resulted, 
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2.71eGo 



hkpVFy/^ — C0S7 \eVmo< 



3/4 



N 
\ n=l 



(9) 



For the fluctuations in the conductance we have de- 
scribed above, terms higher in order than ai^rmn^ were 
not very important. However, when all the N channels 



contributing to the conductance are fully open, the first 
order contribution we have calculated above will be zero. 
Under this condition the second order terms may have a 
noticeable contribution. In this case we can take t = I 
and r = in Eq.|^ which greatly simplifies the derivation. 
It is then quite easy to show that at quantized values 
G = NGo, the contri bution of the sec ond order term in 
ai^r,n„ is CGv oc Ny^ (|a/™„p|ar„„p)- These terms are 
too small to explain the reduction of the depths of the 
minima in the experiment discussed below and will be 
further ignored. 

These higher order terms, however, are not negligi- 
ble when we consider the total conductance of the con- 
tact. The importance of these higher order corrections 
becomes apparent when we compare our theory to the ex- 
periments, and notice that these effects result in a signifi- 
cant correction to the total conductance G. The necessity 
to include these conductance corrections in our analysis 
is a direct consequence of the fact that in the experiment 
we can not measure the bare contact conductance alone, 
as we always measure it in series with the banks. This 
feature in the conductance of quantum po iut con tacts is 
usually referred to as the series resistancetSEacJ. 

The lowest order correction to the average total trans- 
mission probability is given by (Tr[ttt|]) — J2n=i ^"(1 ~ 
Em=i Tm{{\ai,^J^) + {\ar^^ p))). The last term describes 
the path of an electron that is transmitted through the 
contact, scattered back towards the contact in the diffu- 
sive bank, and then transmitted through the contact a 
second time in opposite direction. These processes will 
lead to a smaller conductance than expected for the bare 
contact conductance alone since part of the transmitted 
electrons are scattered back, reducing the net forward 
current flow. 

At higher conductance values, we expect a significant 
contribution of even higher order terms in the return 
probability ai,r„„, to the conductance correction, and 
that hence the lowest order correction used above will 
not suffice. Keeping track of higher order terms, be- 
comes very complicated for many channels. However, 
using random matrix theory an expression for the cor- 
rection to the conductance of a quantum point contact 
connected to diffusive leads has already been derived by 
Bcenakkcr and MelsenEj, 



(G) 
Go 
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l + {g+l)r 



Here g — the reduced conductance of the bare 

contact, where in the theory all channels were assumed 
to be perfectly transmitting. The diffusive scattering in 
the banks is represented by r = Go/Gd, with Gd the dif- 
fusive conductance of the banks. When one makes the 
assumption that the conductance of the banks is large 
compared to the conductance of the contact, Eq. can 
be simplified and rewritten to a form where the correc- 
tion to the conductance effectively becomes a somewhat 
contact dependent series resistance. 



(g + 1)^ 

l + (5 + l)r 



(10) 
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To lowest order Eqs.|l^ and ^ are consistent with the 
correction to the average total transmission probability 
derived from the backscattering above. 



EXPERIMENTAL METHOD 

We have used an MCB to rgftke stable atomic size 
metallic contacts. The techniqueB'H, uses a notched metal 
wire glued onto an electrically insulated phosphor bronze 
bending beam. In the case of sodium a slightly different 
sample fixation method was-*equired due to its high oxi- 
dation rate (see Krans et al&B). The sample is placed in a 
three point bending configuration, which consists of two 
counter supports and a drive at the center of the bending 
beam powered by a piezoelectric element in combination 
with a mechanical screw. By first turning the screw and 
later expanding the piezo by applying a voltage over it, 
we can bend the substrate in a controlled way elongating 
the wire until it finally breaks. The wire is broken at low 
temperature (4.2 K) in an evacuated can to ensure that 
two clean freshly broken surfaces are measured. The volt- 
age applied over the piezoelectric element is linear with 
the elongation of the contact (for further details see Ref 

The conductance measurements are performed by ap- 
plying a 48 kHz, 20 mV amplitude, sinusoidal voltage over 
the contact, which is in series with a 100 fl resistor. The 
first and second harmonic of the voltage over the resis- 
tor are measured, from which we obtain the first (G = 
dl/dV) and second derivative {dG/dV = d^I/dV^) of 
the current with respect to voltage of the contact. The 
conductance is determined with an accuracy better than 
1% for values larger than O.SGq. We use a HP3325b func- 
tion generator to produce the modulation voltage while 
two Stanford Research SR830 lock-in amplifiers at / and 
2/, with a time constant of 10 ms, are used to obtain 
the first and second derivative. 16 bit analog to digital 
and digital to analog converters are used to control and 
measure the piezo voltage. A PC-based controller sweeps 
the piezo voltage up, and while the contact breaks, the 
readings of G and dG/dV are taken through an IEEE 
connection every 100 ms. A full curve of the contact from 
a conductance of over 20 Gq to the transition to vacuum 
tunneling is recorded in about 30 seconds. 

A large number of such curves have been taken for 
gold, silver, copper, sodium, aluminum, niobium and 
iron. In order to avoid anomalously large dG/dV values 
due to unstable contacts near a conductance step, and 
to avoid measuring the average properties of two differ- 
ent plateaus as a result of the finite integration time of 
the lock-in amplifiers (which average G and dG/dV over 
10 ms), only points on a plateau are included through ex- 
clusion of data points, with suitable selection criteria, for 



which the deviation of G and dG/dV with respect to pre- 
vious and consecutively recorded data points is too large. 
After applying this exclusion procedure to each of these 
materials, we have analyzed the results by combining all 
the data and sorting it according to conductance. Then 
a fixed number of consecutive data points were taken 
from which acv = ^{{dG/dVY) - {{dG/dV)^) and the 
average conductance value were determined. With this 
method we obtained aav as a function of conductance, 
in a way which is independent of the number of sampled 
points at a particular conductance value. 



THE NOBLE METALS 
COPPER, SILVER AND GOLD 



The three panels in Fig.|2| show the measurement of 
the differential conductance, obtained with a small (< 
0.35 mV) modulation voltage, against bias voltage for 
three gold atomic size contacts. In each case two curves 
are plotted, one for increasing and one for decreasing bias 
voltage, showing the reproducibility of the behavior. The 
bias voltage over the contact determines the energy of the 
electrons injected into the banks and hence modifies the 
electron interference resulting from electrons scattered 
back towards the contact in the banks. This change 
in the interference gives rise to the fluctuations shown 
in Fig.^. In our experiments described below the volt- 
age dependence of the conductance is determined with a 
modulation amplitude of 20 mV, i.e., the average slope 
of curves such as those presented in Fig.^ is determined 
over a bias voltage range of ±20 mV. Note the small am- 
plitude of the fluctuations in Fig.|^b. We will argue later 
that this is an example of the reduction of acv for the 
conductance of gold contacts with value near 1 Gq due to 
the v/T2(l - r„) factor in Eq.|. 
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FIG. 2. Plotted in the three panels is the differential con- 
ductance dl/dV as a function of bias voltage, measured 
with a modulation amplitude < 0.35 mV, for three differ- 
ent gold contacts with G (a) ~ O.88G0, (b) ~ 1.02Go, and (c) 
~ 1.65Go. For all three curves the vertical-scale spans 0.12 
Go. 
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An example of the typical conductance and dG/dV be- 
havior when breaking a gold contact for a constant mod- 
ulation amplitude of 20 mV and zero bias is presented 
in Fig.^. The steps and plateaus in the conductance 
correspond with atomic rearrangements and elastic de- 
formation resgectively, as the contact is pulled apart and 
finally breakg3. At each step in the conductance we find 
a corresponding step in dG/dV. Even tiny steps in the 
conductance, such as between 7 and 8 Go, can produce 
dramatic jumps in dG/dV. Changes in electron path 
lengths of the order of the Fermi wavelength (which is 
the atomic scale for metals) occur at these steps in the 
conductance as a result of atomic rearrangements, and 
randomly change the resulting electron interference. The 
continuous change of dG/dV during elastic deformation 
of the contact along a plateau results from the gradual 
elongation of the electron path lengths and hence in a 
gradual change of the resulting interference. 
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FIG. 3. Typical example of the simultaneous measurement 
of voltage dependence of the conductance dG/dV and the 
conductance G, as a function of piezo voltage Vp for gold 
measured with a constant modulation amplitude of 20 mV. 
The graph includes vertical gray lines which show that the 
steps in both quantities coincide. Two plateaus have been 
enlarged and offset to show the tiny steps in the conductance. 
The open squares represent the points excluded from the en- 
semble average by the selection procedure. The elongation of 
the contact is linear with Vp and 10 V corresponds to about 
1 nm. 



In Fig. H the open squares represent the points at steps 
in the conductance and dG/dV which have been excluded 



from the statistical analysis by the selection procedure 
discussed above. As can be seen in the figure, the ex- 
cluded data consists exclusively of the last and first points 
on a plateau, together with points which lie between two 
plateaus as a result of the finite integration time of the 
lock- in amplifiers. 

Fig.^ shows the distribution of values measured for 
dG/dV in a particular range of conductance collected 
from 3500 individual curves similar to the one presented 
in Fig. ^. The distribution at 1 Go is clearly much nar- 
rower than the other two at non-integer values. This 
is statistical evidence for what was already observable in 
measurements of the bias dependence of the conductance 
for a single contact in Fig.gb. We will argue that the 
narrow distribution can be explained by a suppression 
of the conductance fiuctuations (Eq.^), as a result of Ti 
being approximately equal to 1 and all other T„ w 0. In 
the presentation of the data below we concentrate on the 
width of these types of distributions, ggVi determined 
for a fixed number of data points. 
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FIG. 4. The distribution of dG/dV values in a particular 
conductance range collected from 3500 individual curves for 
gold such as the one presented in Fig. 3. The conductance 
range of the three curves roughly corresponds with the con- 
ductance of the dl/dV curves in Fig 2. (a) G = 0.9 ± 0.05Go 
(b) G = 1.0 ± O.OSGo (c) 1.6 ± 0.05Go. The dotted curve in 
(a) represents a Gaussian fit of the data. 

The shape of the distribution curve shown in Fig.^ 
is much sharper around the dG/dV = value than a 
Gaussian distribution. The tails of these curves also 
deviate from Gaussian behavior. These deviating fea- 
tures are analogous to the peaks calculated and mea- 
sured in the distribution of parametric derivatives (e.g. 
thermopower)[€>f quantum dots with single mode ballistic 
point contactg23. The origin of this cusp at zero ampli- 
tude is the limitation of the range over which the dif- 
ferentiated parameter can vary in value. At both of its 
maximal values the parametric derivative is zero, leading 
to an enhancement of the statistics at zero amplitude. 

In the upper panels of Fig.||, we present the mea- 
sured acv for the noble metals copper, silver and gold. 
The data points for ~ 3000 individual curves such as the 
one in Fig.|^ were sorted as a function of conductance. 
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From this total collection of data points the root mean 
square of dG/dV was calculated for groups of 300, 2000 
or 2500 successive data points, depending on the density 
of points available. This total collection of data points 
was also used to calculate the corresponding conductance 
histograms plotted in the lower panels of the figure. 



Copper Silver Gold 




G (2e'/h) G (2e'/h) G (2e'/h) 

FIG. 5. (JGv (top) and conductance liislogram (bottom) 
against G for (a) copper 3000 curves (b) silver 2400 curves 
(c) gold 3500 curves. The averages have been obtained from 
300 data points for each of the open circles and 2000, 2000 
and 2500 data points for each of the solid squares for copper, 
silver and gold respectively. The solid curve in each case show 
the behavior of a single channel opening at a time. The ver- 
tical gray lines represent the corrected integer conductance 
values. 

The electronic properties of these three noble metals 
are very similar, which is reflected in the similar behavior 
we obtain for ugv as a function of conductance. Min- 
ima in ugv near 1, 2 and 3 Go can be observed in all 
three cases. The minima, however, are most pronounced 
for gold which even has a small dip near 4 Go- Another 
important similarity, as is apparent from the peaks in 
the conductance histogram, is the preferred values for 
the conductance just below I, 2 and 3 Go for all three 
materialsll3T3. 

When comparing the experimental results with our 
model it is important to note that a given value for 
G = Go X]^^=i '^S'li be constructed in many ways from 
a choice of transmission values {T„}. The experimental 
values for (tgv are, therefore, an average over impurity 
configurations and transmission values. Assuming the 
averages are independent, we can compare the data with 
various choices for the distribution of the transmissions. 
In each upper panel of Fig.|^ the solid curve depicts the 
behavior of Eq.|| for a single channel opening at a time, 
i.e. in the interval G/Gq € {0, 1} there is a single chan- 
nel contributing to the conductance with G = GqTi, in 
the interval {1,2} there are two channels with one fully 
open G = Go(l + 12), etc.. The curves in the figure have 
been scaled so they best fit the data. From this scal- 



ing an estimate for the mean free path can be obtained 
when a reasonable value rangec3 for the opening angle 
7 = 35° — 50° is assumed. For copper, silver and gold, 
we obtain a value of = 3 ± 1 nm, Ig = 4 ± 1 nm and 
If, — 4:± 1 nm respectively. 

For gold the description of the experimental data with 
a single channel opening at a time works surprisingly 
well In particular, for the minimum near IGo, and for 
the fact that the maximal values between the integer 
conductance values are all nearly equal. The main dis- 
crepancy is that the minima become less pronounced for 
higher conductances. The well-developed structure ob- 
served in acv, with a dependence which closely follows 
the ^n(l ~ ^n) behavior of Eq.^, demonstrates a 
property of the contacts which vae refer to as the satura- 
tion of the channel transmissionB: There is a strong ten- 
dency for the channels contributing to the conductance 
of atomic-size contacts of gold to be fully transmitting, 
with the exception of one, which then carries the remain- 
ing fractional conductance. 

For copper and silver the amplitude of the data in- 
creases together with the degradation of the minima. 
These two metals also exhibit the saturation of the chan- 
nel transmission effect, but clearly not as rigorous as for 
gold. This reflects itself in the estimates we can make for 
the contribution of an additional channel at the first three 
minima. Neglecting the small contribution of the higher 
order terms in a;^r,„„, these are for copper 2%, 12% and 
15%, for silver l%ril% and 18% and for gold 0.5%, 6% 
and 10%. The concept of the saturation of transmission, 
channels is consistent with the model of Cuevas et al.U 
and other recent experimental work, which shows that, 
for gold, the conductance f^r-G = 1 Go of a single atom 
is carried by a single modeElO. 

The minima in (jgv lie at values for G below the inte- 
ger conductance values. This shift is due to the scatter- 
ing of transmitted electrons back to the contact, which 
apart from a fluctuating first order contribution in 
which determines dG/dV, also gives rise to a shift in (G) 
when contributions to second order in ai,r^„ are taken 
into account. Ideally we would like to plot acv as a 
function of with T„ the transmission probabil- 

ity of mode n of the bare contact, but the bare contact 
is always measured in series with the diffusive banks. In 
order to correct for the backscattering to lowest order, 
the theoretical curves have been plotted as a function of 
G = GoEtiTnil - E™=ir™((|a,_p) + (|a._n)). 
The (|a;,r„„ P) have been adjusted for optimal agreement 
with the experimental minima. The vertical gray lines in 
each figure represent the corrected integer conductance 
values using this lowest order procedure. For gold and 
silver this value is comparable (|az,r,„„P) = 0.005 ver- 
sus 0.004 respectively, and hence is consistent with the 
similar amplitude of acv observed for both materials. 
Copper has a somewhat larger amplitude for acv , in ac- 
cordance with a larger shift in the conductance minima 
with (|a/,r„„P) — 0.014. We can relate this value for 
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the conductance correction to a rough estimate for the 
mean free path by equating (|a/.r„,„P) to the integral of 
Eq.|^ over all path lengths. As lower integration limit 
we have taken the typical shortest path time Tg, for the 
upper limit infinity and for the opening angle we have 
assumed the typical range 35° — 50°. The value for 1^ 
we obtain using this method is 4 ± 1 nm, 7 ± 1 nm and 
6 ± 1 nm for copper, silver and gold respectively. This 
is quite close to the mean free path derived above from 
the amplitude of acv and therefore in accordance with 
our model. The values for the mean free path we obtain 
are much shorter than what is normally found for bulk 
samples, and can probably be attributed to surface scat- 
tering near the contact. Assuming surface scattering is 
indeed responsible, an important property of the mean 
free path which we neglect here is that Zg will not be 
a constant as a function of the conductance, but rather 
increase as the contact diameter becomes larger. This 
size dependence of the mean free path is not expected to 
be very significant in the range of validity of our model, 
where we assume a contact diameter d < L < 

The relatively short 1^ we obtain is responsible for 
the correction to the quantized conductances, but it is 
too long to hold backscatting responsible for the mea- 
surement of the significant frequency with which non- 
quantized values are measured. Also, if scattering is held 
primarily responsible for reducing the conductance from 
for instance a perfect conductance of 2 Go to 1.5 Go, then 
it is not unreasonable to assume that contacts with a 
perfect conductance of 1 Go are reduced to 0.5 Gq with 
a probability of the same order of magnitude. This is 
not observed experimentally at low temperatures, as con- 
tacts with a conductance of 0.5 Go occur more than 500 
times less frequently for silver and copper than contacts 
with a-conductance of 1.5 Go (The formation of atomic 
chain^Sl reduces this ratio to about 20 times in the case 
of gold, since the conductance of the chains is quite sen- 
sitive to distortions making contacts with a conductance 
of 0.5 Go occur with an enhanced frequency). If, on the 
other hand, one assumes that contribution from tunnel- 
ing, due to for instance geometrical considerations, are 
more important, the appearance of non-quantized val- 
ues above 1 Gq finds a natural explanation. The for- 
mation of geometries with a conductance smaller than 
1 Go is highly unlikely since the smallest contact geom- 
etry is that of a single atom with conductance 1 Go and 
when the contact breaks, the banks relax back prevent- 
ing high transmission probability tunneling contributions 
from contributing. The latter process is usually referred 
to as the jump to tunneling^. 

An important feature for all three noble metals is that 
the minima in acv near 2 Go do not coincide with the 
respective peaks in the histograms. The minima lie at 
the expected conductances based on the backscattering 
amplitude we require to consistently fit all the minima, 
the so called corrected integer conductance values. The 
second peak in the histograms clearly are located at lower 
conductance values. We propose that this discrepancy is 



caused by favorable atomic configurations which have a 
bare conductance smaller than 2 Go, and thus give rise to 
a peak in the histogram below the corrected conductance 
value for 2 Go . 



TEST OF THE MODEL 

Apart from the set of transmission values {Tn} the 
equation contains two free parameters, le and 7, over 
which we do not have any experimental control. The de- 
pendence on the modulation voltage Vmod, however, is a 
parameter which we can verify, assuming that all other 
relevant parameters are independent of the applied volt- 
age. This is a reasonable assumption for the modulation 
amplitudes at which we measure, as parameters like {T„} 
vary on the scale of Ep, and and 7 are not expected 
to change on the scale of lOOmV. Assuming further that 
eVmod > Tt-It^ then the product crcyV^^^^ should be con- 
stant for all Vraod- 

In Fig. |[ c^GvV^^^ has been plotted against G, for 
Vmod = ^^^ 20, 40 and 80 mV. Within the experimental 
accuracy no modulation voltage dependence is observed, 
as all four data sets coincide very well. This would not 
be the case unless the power of the modulation ampli- 
tude dependence is close to 3/4. Using a procedure that 
calculates the minimal difference between the six com- 
binations of experimental curves, which have been mul- 
tiplied by their respective modulation amplitudes to a 
power which is the free parameter, we find this power to 
be 0.71 ± 0.06 in good agreement with the 3/4 predicted 
by the theory. 




0.00 -F ■ 1 ■ 1 ■ 1 ■ 1 1 

12 3 4 

G (2e'/h) 

FIG. 6. acv multiplied by the modulation amplitude to 
the power 3/4 against conductance G for gold measured at 
(squares) 10 mV, (circles) 20 mV, (up triangles) 40 mV, (down 
triangles) 80 mV. 

The mechanism used to describe the fluctuations in the 
conductance above also produces fluctuations in other 
transport properties, notably the thermopower. Measur- 
ing the thermopower of atomic size metallic contacts re- 
quires a completely different experimental method and is 
performed on an energy scale much smaller than that nec- 
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essary for determining acv ■ The experimental resultfl, 
however, have been successfully described by a theory 
based on the same principles as those presented above. 
The predicted theoretical relationship between the stan- 
dard deviation of the thermopower as and acv is given 
by, 



THE ALKALI METAL SODIUM 



2.71e2GoELi^n 

ckB{kBeY/\eV,nodf'^' 



(12) 



where the numerical constant c = 5.94. The only 
parameters in the above relation are, 9, the tempera- 
ture at which the thermopower measurements were per- 
formed, Vmod, the modulation amplitude for the conduc- 
tance fluctuation measurements and the total conduc- 
tance GoX]n=i-^"- unknown parameters in Eq.j^, 
the set of transmission probabilities {Tn}, Ze and 7 can- 
cel out, and the scaling relation provides an independent 
test of the experimental data. 
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FIG. 7. Comparison of the standard deviation of the 
thermopower as and the standard deviation of the volt- 
age dependence of the conductance ogv by plotting (Q) 
2.7le'^(Js/ckB{kBeY/^{eVmodf^'^ and (□) ggv/G as a func- 
tion of conductance. 



By comparing acv to measurepients of the ther- 
mopower on atomic size contactsE2l using Eq.|l^ we 
can effectively extend the energy range over which 
we test our model to an order of magnitude smaller, 
and test the experimental procedure against a com- 
pletely independent method. In Fig.|^ we have plot- 
ted 2.ne^as / ckB{kBOY'^{eVmodf'^ and ugv/G as a 
function of conductance. We have set 6* = 12K and 
Vmod = 20 mV in accordance with the experimental con- 
ditions. Excellent agreement is obtained between both 
experimental methods without any free parameters. We 
interpret this as a successful test for the validity of the 
principle on which our theoretical analysis is based. The 
dependence on the opening angle 7, however we regret- 
fully cannot verify experimentally as we have no control 
over the contact geometry. 



Sodium also is a monovalent metal but its histogram 
determined from 1800 curves (lower panel Fig.||), is com- 
pletely different from that observed for copper, silver or 
gold. The statistically preferred conductance values are 
observed as peaks in the histogranxj near 1, 3, 5 and 
6 Go rather than near 1, 2, SGq. This series of peaks in 
the histogram at 1, 3, 5 and 6 Go have been interpreted 
as resulting from the quantization of the conductance in 
a cylindrically shaped nanowire. The histogram peaks 
are very sharp, and in the 1800 curves measured almost 
no data is obtained between and 1 Go and between 1 
and 2 Go. For this reason no points for cjgvi determined 
from the same 1800 curves, are presented in these ranges 
(upper panel Fig. Even with these points absent, the 
ugv measured for sodium is distinctly different from that 
observed for the noble metals. In acv we observe def- 
inite minima near 3 and 6 Go and although there has 
been no data measured in the ranges < G < 1 Go and 
1 < G < 2Go, the value of acv at 1 Go is small making 
it a very probable location of a minimum. Since there 
is no data below 2 Go, we cannot exclude that there is a 
small minimum at 2 Go . 

The histogram peaks coincide with the minima in acv , 
with the exception of the peak near 5 Go. The ab- 
sence of a minimum at 5 Go is at first surprising. When 
one considers that in a conductance histogram for a 
model of a three dimensional cylindrical contact based 
purely on a free electron gas, the peak at this value is 
found to be iiearly absent due to smearing by tunneling 
contributionsEj, it is striking that a histogram peak is 
there at all. Also, unlike the other peaks, the one below 
5 Go and also the small one above 2 Go do not coincide 
with the corrected integer conductance values. We pro- 
pose that these two peaks result from favorable atomic 
configurations, which are sampled more often than other 
conductance values while making a histogram, but do 
not result from stable quantized conductance values de- 
termined by an integer number of nearly open channels. 

As in the case of the noble metals there is a system- 
atic shift of the position of the minima in ugv to lower 
conductance values. The corrected integer conductance 
multiples are shown in Fig.|^ as vertical gray and black 
lines depending on whether they coincide with the min- 
ima in ugv or not. For the correction to the integer 
conductance values in Fig.|| we have used Eq.|l^ with 
r = 0.015. This value corresponds to a series resistance 
of 200 Q, for g ^ 00 and 400 for (7 = 1 as is evident from 
Eq. 
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Eq. |10| is used, rather than the first order correc- 
tion applied to the noble metals, because the conductance 
extends to larger values where higher order terms in the 
backscattering amplitude become relevant. 
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Sodium 




G (2e'/h) 

FIG. 8. cjgv (top) and conductance histogram (bottom) 
against G for 1800 sodium curves with each sohd square repre- 
senting the statistics on 1000 data points. The vertical black 
and gray hnes indicate the corrected integer conductance val- 
ues for which the histogram peaks respectively do and do 
not correspond with minima in aav- The curves depict the 
behavior of a hyperbolic constriction in a three dimensional 
electron gas with circular aperture, with ( — ) opening angle 
7 = 60° and mean free path le = 4.4 nm, (■ ■ ■) 7 — 45° and 
le = 5 nm. 

We want to compare these results with the simplest 
possible model, which neglects the atomic character of 
the contact and only takes the cylindrical symmetry and 
finite length of the contact into account. For this purpose 
we have made use of the calculations by Torres et alS3. 
The model consists of a free electron gas confined by hard 
wall boundaries, which have the form of a hyperboloid. 
The differential equations for this system can be solved 
numerically, from which the transmission probabilities r„ 
for each mode as a function of the contact diameter can 
be obtained. This makes a direct evaluation of Eq.^ pos- 
sible, where the only remaining adjustable parameter is 
the mean free path. The opening angle which describes 
the shape of the hyperboloid and thus the dependence of 
the mode transmissions on the contact diameter, is the 
same as the one which enters in Eq.||. We have added two 
acv curves calculated for such a system to the graph. For 
curves in the opening angle range from 60°, with mean 
free path Zg — 4.4nm, to 45°, with = 5nm, we find 
reasonable agreement between various ranges in the data 
and the theoretical curve. 



The differences between the calculated curve and the 
measured data can be attributed to the averaging over 
many contact geometries and thus over a range of 7 val- 
ues. Also, the smearing in the conductance (G) due to 
the ensemble average of defect configurations is not in- 
cluded. This property will make the minima less deep 
and sharp but will hardly influence the maxima. 

Another feature of the calculated curve, which can also 
be recognized in the 

measured data, is that the minimum in the experimen- 
tal and calculated (tqy below 6 Go does not coincide ex- 
actly with the corrected quantized value for the conduc- 
tance, even after application of the same series resistance 
we have used to compensate the shift in the histogram 
peaks. In other words, when we ignore the series resis- 
tance correction, the model predicts the minima to be 
shifted above the integer values. This is a direct result 
of the signiflcant tunneling contributions for the opening 
angles we have used to model the contact when just open- 
ing a mode in combination with the asymmetry of the 
dependence of ctgv on the mode transmission. This sys- 
tematic shift becomes more pronounced with larger open- 
ing angles, larger conductance values and the presence of 
degenerate modes. The value of the opening angle we ob- 
tain for sodium from the theoretical curves, 7 = 45° — 60° 
is comparable but somewhat larger than the typical es- 
timates KLade for the opening angles of atomic size gold 
contactsEa. 



THE TRIVALENT METAL ALUMINUM 



The statistically preferred conductance values for alu- 
minum are shown in the lower panel of Fig.^. The clear 
peaks, evident in the histogram below 1 and 2 Gq and a 
weak bump above 3 Go are in accordance with previous 
measurementaH3. The peaks are less pronounced, but at 
first glance similar to those observed for gold, silver and 
copper. The most important discrepancy between the 
monovalent metals and aluminum is that in the latter 
case the first peak is broader and clearly displaced below 
IGo. 

The measured crcv for aluminum, presented in the up- 
per panel of Fig.^, is completely different from the be- 
havior observed for the noble metals copper, silver and 
gold. The clear minima at 1, 2 and 3 Go have been re- 
placed by a slight dip at IGo. In order to understand 
these measured features it is important to realize that a 
single aluminum atom has a conductance close to Go but 
admits three conductance channelsElEEl. It is thus not sur- 
prising that the behavior associated with the saturation 
of a single partially open channel is not observed for alu- 
minum. The histogram peaks observed for this trivalent 
material can thus be attributed to another mechanism. A 
likely candidate is favorable atomic configurations, which 
are probed more frequently than others. 
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Aluminium 




G (2e 7h) 

FIG. 9. (JGv (top) and conductance histogram (bottom) 
against G for 2800 aluminum curves with each sohd square 
representing the statistics on 2000 data points. The curves 
in the graph have been labeled with a number and represent 
the contribution to <jgv of (1) a single channel opening, a 
random distribution over (2) two, (3) three, (4) four, (5) five 
and (6) six channels (see text). 

In Fig. ^ a series of curves are included which show the 
behavior for a single channel, and a random distribution 
of two, three, four, five and six channels. The curves 
have been generated by calculating the square root of 

O"ma;r/o -/o ^C^l ' , ■ • , TV ) En=l (1 " ) dTi . . dT^r 

where P is the probability distribution giving an equal 
probability to every transmission value for each T„ un- 
der the constraint that ^« — G/Gq. These curves 
have been labeled 1, 2, 3, 4, 5, and 6 respectively, and 
have all been scaled with the same amplitude. The dips 
in the calculated curves with a random distribution of 
two or more channels results from the property that, in 
a random distribution, at multiples of the conductance 
quantum, there is a finite probability to encounter some 
r„ — (0 or 1), for which their contribution to aov — 0- 
This effect however, and hence the dip, becomes less pro- 
nounced with an increasing number of channels. For the 
lowest conductances (< 0.5 Go) one can observe that the 
behavior most closely follows that of a single channel. 
This is expected as for small conductances in the tunnel- 
ing regime a single channel is expected to dominate the 
conductance and, hence, the behavior in acv- As the 
conductance increases the behavior becomes more like a 
random distribution over an increasing number of chan- 



nels. The random distribution of channels which we have 
introduced serves only to illustrate that the small dip at 
1 Go appears for a limited number of channels even with- 
out the saturation of channel transmission effect, and 
that the gradual increase in acv is a direct result of 
more channels contributing for larger conductances. The 
actual behavior of the transmission channels is proba- 
bly not completely random as can be judged from the 
theoryQ which shows that there is usually one dominant 
channel and two smaller ones for a single aluminum atom. 
Nonetheless the curves reproduce the evolution of acv as 
a function of conductance with reasonable accuracy. 



THE TRANSITION METALS 
NIOBIUM AND IRON 



The measured cjqv (upper panel) and histogram (lower 
panel) for 2400 niobium curves, recorded at a tempera- 
ture of 10 K in order to avoid effects of the supercon- 
ductivity on the voltage dependence of the conductance, 
have been plotted in Fig.jlCja. For iron the measurements 
of acv (upper panel) and the histogram (lower panel) 
recorded for 700 curves are presented in Fig.nT 



Niobium 



Iron 




G (2e7h) G (2e'/h) 

FIG. 10. (a) GOV (top) and conductance histogram (bot- 
tom) against G for 2400 niobium curves measured at 10 K 
with each solid square representing the statistics on 1000 data 
points, (b) acv (top) and conductance histogram (bottom) 
against G for 700 iron curves with each solid square repre- 
senting the statistics on 500 data points. 

Both sc?-metals, niobium and iron, show completely 
different features when compared to the other materi- 
als we have discussed so far. When compared with each 
other, the measurements for niobium and iron are so sim- 
ilar that they nearly are indistinguishable. For both met- 
als, (Jgv increases strongly from to 1 Go and above this 
conductance value, increases only slightly. The dip ob- 
served in (jQv for aluminum is absent and the increase in 
UGv with conductance is much smaller than was the case 
with aluminum. From the single peak in the histogram 
we can deduce that both materials have a statistically 
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preferred conductance value just above 2 Go- This peak 
is expected to be the result of the reproducibility in the 
conductance of the last plateau consisting of a single nio- 
bium or iron atom. In the case of niobium this value is 
in excellent agreement with the measured and calculated 
conductance value for a single atom of niobium (2 to 
3 GoQu) . For iron these calculations and measurements 
have not yet been performed. 

The completely random distribution used to describe 
the behavior of aluminum clearly cannot reproduce the 
measured aov for niobium and iron. If one considers 
distributions for the transmission probabilities that are 
closer to the five calculated transmission probabilities 
contributing to a single-atom contactQ, the experimen- 
tal behavior up to 3 Go can in principle be simulated. 
The number of free parameters in such an analysis, how- 
ever, makes such an ad hoc procedure quite meaningless. 
More theoretical work must be performed to provide an 
approximation for the range of transmission channel dis- 
tributions that should be considered in order for the mea- 
surements to be reliably related to a general trend. The 
fact that the measured acv above 1 Go is almost not 
dependent on the total conductance suggests that the 
number of partially open channels contributing to the 
conductance does not depend on the contact size. 

The histogram we have recorded for iron is different 
from the histogram constructed from 80 iron contacts 
recorded under—ambient conditions in the magnetically 
saturated stateEj. Possibly the different temperature at 
which the experiments are performed or the magnetiza- 
tion state are responsible for the discrepancy. The influ- 
ence of these experimental conditions should be studied 
in more detail in the future. 



CONCLUSIONS 

With the technique of ensemble averaged conductance 
fluctuations we are able to measure the average prop- 
erties of the conductance mode evolution of atomic size 
contacts. We have successfully tested the modulation 
voltage dependence of the theory, and can relate acv 
to the standard deviation in the thermopower without 
any free parameters. An important property of the mea- 
surement of the conductance fluctuations is that it is not 
dependent on preferential contact configurations. 

Sodium, the most free electron like material studied, 
exhibits electronic behavior which can be reproduced rea- 
sonably well with a hyperbolic constriction in a three 
dimensional electron gas with a circular orifice. The con- 
ductance histogram of this material however, contains 
contributions from such a circular constriction in a three 
dimensional electron gas as well as other peaks possibly 
resulting from favorable atomic configurations. The con- 
ductance properties of the other monovalent materials 
such as gold, silver and copper seem to be best described 
by the tendency of the conductance channels to open one 



by one, a property which has been called the saturation 
of channel transmissioro. The conductance histogram 
of these materials contains features which seem to coin- 
cide with the evolution of the conductance modes, but 
particularly the second peak in the histograms are also 
determined by other statistical (probably atomic) prop- 
erties of the contact. For aluminum, niobium and iron we 
find the behaisiopior ugv we expect based on the atomic 
orbital modelElu'tJ. The conductance histograms of these 
three materials seem to be dominated by the statistical 
distribution of atomic configurations. 

The concept of saturation of channel transmission has 
been introduced in Ref. ^, to make a marked distinction 
between the properties of the conductance modes we ob- 
serve here, and the statistically preferred conductance 
values which in the literature are generally referred to 
as conductance quantization. Indeed, both are based on 
the quantum mechanical Landauer-Biittiker formalism, 
but the latter assumes there is a statistical preference for 
contacts with quantized values as a result of this formal- 
ism. Since it is not clear to what extent favorable atomic 
configurations are responsible for the histogram peaks, 
we feel a sharp distinction should be made between re- 
sults that can be influenced by favorable atomic positions 
possibly mimicking the features of conductance quantiza- 
tion, and results which truly probe the electronic prop- 
erties of the contact. With the saturation of the channel 
transmission we wish only to describe the evolution of 
these modes, but we do not rule out that conductance 
quantization may prove to be anj-important factor which 
influences the contact formationtil. However, in view of 
the results presented in Figs.^ and |l^ together with the 
arguments presented above for the other materials, we 
should be aware that peaks in a histogram by themselves 
give no unambiguous information about the actual com- 
position of the conductance modes. 

For the correction to the bare contact conductance due 
to scattering near the contact, we find that this correc- 
tion is correlated to the amplitude of the conductance 
fiuctuations and hence the elastic mean free path. This 
provides strong experimental evidence that these types of 
scattering effects are indeed responsible for the so-called 
series resistance. 
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APPENDIX 

Incorporating a finite modulation voltage into the ex- 
pression for (7gv is somewhat more technical. In analogy 
with Eq.| we can write for the energy dependent part of 
the current, 



eV N 



rev J^" 

/ Vr„2Re[r„a,„„(-£;) 
Jo , 



+r'r,ar^SeV - E)]dE. 



(13) 



We must include the modulation voltage F = Vq + 
KnorfSin(a;oT') in the argument for the backscattering 
amplitude, which we write as an integral over contri- 
butions of path traversal times r. Hence ai,r„„ {E) = 
J ai^r„„{'r)exp[~iET/'h]dT where we assume the domi- 
nant energy dependence is the phase factor. We will first 
consider the contribution to the current of a single path 
with time t: SI{t). At a later stage we will perform the 
integration over different paths. Evaluating the integral 
over E in Eq. O gives, 



N 



51 (t) = const. + ^2 V r„- 

h ^ — ^ 7- 



(14) 



n=l 



xRe[-ir„a,„„(T)e*'=^«^/''e+*'=^"<-'™"('^''^')/'' 
+«^>r„„ (T)e-*'=^"^/''e-*'=^'"°''"'"('^°^')/'*] . 
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We are measuring the ac component of the current at 
twice the modulation frequency, this is equivalent to mea- 
suring the second derivative of the current with respect 
to voltage. Therefore we expand the exponential term 
Eq. 14 in harmonics of the modulation frequency, 



n— — C30 



Jn{eVraodT/h) 



= Jo{eVmodT/h) + 2ism{uJoT')Ji{eVjnodT/h,) 

+ 2cOs{2LUoT')J2{eV,nodT/h) + (15) 

where Jn{z) is the n'th Bcssel function. In the last step 
we used J_„(— z) = (— 1)" J„(z). We are particularly in- 
terested in the last term in Eq.|l^, which we will use to 
obtain the part of the current which is proportional to 
twice the modulation frequency (cos(2ci;o''"'))i 

5I2ujo{t) = --^4^T„-cos(2woT-')^2(eKnodT/?i) 

n— 1 

xRe[-irnai^ST)e''^°^'^ + ir[,ar„^{T)e-"'^°^'\ (16) 

In general, the Taylor expansion of the current for small 
modulation amplitude is 

/ {Vo -f VraodSH^t)) = /(Vo) + (^^^ Sin{ujt) 

1 / 9^/ \ 

"4 Vmn) '^^odC0s{2ujQT') + .... (17) 
Combining Eqs.|l|, |l| and |l| we obtain, 



4 2e ^ fi 

TT2 7-2 ^ Tn-J2{eV,nodT/h) 

^mod " n=l ^ 



X (-*^«ai„„(T)e*^^''"/'' + ir>,„„(T)e-*^^°"/'^ + 
«„Me""'^°^/'' - ^<*<_(T)e+^'=^°^/'')^^T. 



We now can continue as before and calculate the corre- 
lation function. We assume only terms with the form 
^'nn (''")'^/* (''') or ar„„ (r)a*^^ (r) contribute to this func- 
tion as all other combinations of reflection coefficients 
and combinations of different path times will average to 
zero, and introduce the brackets which represent averag- 
ing over different impurity configurations. 



^GV 



16 



\/4 



1- 



N 

E 



; T„2(l - T„) (^^) Ji{eV„.odT/h) 



X (■^)«r„„ i'^) + (•r)a*„„ (T))dT. 

Both sides of the contact have the same average 
properties, and we can write for {ai^^{T)al (r)) = 
(ar„„(T)a*^^^(T)) ~ Pci{t). Filling in the expression for 
the classical return probability (Eq.^) and substituting 
X = eVmodr/h gives us. 



^GV 



^morf \ / n=l 



2V3^fc2,i:>3/2(l-cos7 



T/ \ 



Numerical integration of the integral yields the value 
0.03385 and filling in D = v'^Te/S finally results in Eq.||. 
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